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ABSTRACT 



A method of deriving a probability distribution function 
to approximate the distribution of times required for the 
conduct of military operations is developed, A procedure for 
collecting the necessary data is suggested, and a detailed 
description of an algorithm for manual generation of random 
times is given. Thus the capability of introducing random 
times to simulate time delays in manual war games is pro- 
vided. 
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PREFACE 



War games are currently used to a wide extent to simu- 
late the interactions of opposing forces in the conflict 
situation* The validity of conclusions based upon the re- 
sults of such simulation depends directly upon the quality 
of the inputs to the games. One important input is the 
time required by forces to accomplish certain missions or 
conduct specific operations. Many games presently utilize 
a method of fixed time delays which remain constant for 
specific operations throughout the play of the game. De- 
pending upon the objectives of the war gaming effort, this 
may well be adequate. Under certain conditions, however, 
it is desirable, if not mandatory, to utilize some method 
of random time delays. 

Where an attempt is made to follow the latter course, 
the current trend is to introduce times generated from 
either a Uniform Probability Distribution or a Triangular 
Probability Distribution. The use of these distributions 
is quite arbitrary, based upon ease of computation rather 
than degree of approximation to reality. In effect, they 
are used only to inject an element of randomness rather than 
to simulate the random manner in which real times vary. 

In introducing random time delays to a manual war game, 
there are two basic tasks to be accomplished. One must first 
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derive some probability distribution which varies in a 
manner similar to the manner in which the real times being 
simulated vary. Having done so, one must then devise some 
method to generate random numbers from this distribution, 
which will be used to simulate random times. 

This thesis attempts to present a method of accomplishing 
these two tesks which offers an advantage over the methods 
presently in use. A procedure is suggested, where in the 
absence of recorded or experimental results, professional 
estimates of the required times may be used to determine an 
approximating probability distribution. A method of generating 
random times from this distribution is presented which repre- 
sents a compromise between ease of computation and degree of 
simulation. 

Chapter I considers the overall problem of determining 
the specific probability distribution which best describes 
a real life event, and describes the difficulties inherent 
in solving this problem. 

Chapter II describes a method of collecting certain 
minimal data about the distribution of times required for 
the conduct of military operations. It further demonstrates 
how this data may be used to develop a distribution function 
which approximates the distribution of the real times. 

Chapter III presents an algorithm for generating random 
time delays based upon the approximation resulting from 
Chapter II. 
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INTRODUCTION 



SECTION 1 

War gaming is an attempt to build a workable and 
meaningful model of the interactions between opposing forces 
in a conflict situation, and in this paper should be dis- 
tinguished from war exercises in which troops are actually 
deployed in the field. Whether on a map, computer, or 
playing board, war gaming is an abstraction or idealization 
of reality and is in no sense an attempt to duplicate reality, 
which attempt would require the philosopher’s "perfect” 
knowledge of the universe. The level of this abstraction 
or idealization has varied throughout the history of war 
gaming from the highly idealized games of war chess (circa 
I 6 I 4 .Z 4 . ) through Prussian rigid Kriegsspiel (circa 1865) with 
its complicated and elaborate rules and tables and the various 
forms of free Kriegsspiel (circa 1875) with their heavy 
reliance on control and director groups, free play and flexi- 
bility. ( 1 ) 

Historically, rigid Kriegsspiel was an attempt to 
move the war game out of the parlor and closer to real oper- 
ations by using maps, actual rates of fire, actual movement 
times and a set of rules defining how assessments were to 
be made. In its later development, however, it suffered 
from this same elaborate structure of rule books, tables and 
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formulas necessary for the assessment of play, so that play 
lost a great deal of flexibility and consumed considerable 
time. As a reaction against this complication, free 
Kriegsspiel was developed in which the assessments of play 
and guidance for the game were made by a director or control 
group staff, materially reducing the necessity for elaborate 
formal rules but placing heavy burdens on the operation and 
method of selection of the control group. 

Modern war gaming preserves this distinction between 
’’rigid" and "free" war games quite well. Until recently, 
most war gaming was of the "free" play type mainly because 
of its relative ease of play, its flexibility and the ability 
of the control group to determine the depth of detail to be 
investigated. With the development of high-speed digital 
computers to handle the extensive tables and rules, however, 
"rigid" type games have again come into their own. So the 
modem manual war game (say, the CPX) is a direct descendant 
from rigid Kreigs spiel. 

As there is no such thing as the "right" type of war 
game but only the 'toost useful" for a particular purpose, 
these two extreme forms of the war game are often combined 
in varying degrees and as a group are called composite war 
games. The games which are more formal and rigid are par- 
ticularly well suited for investigations of operations in 
which a great deal of detailed information is available and 
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the interrelationships of the various components of the 
game are well understood. These games can often provide 
useful information for current planning, research and evalu- 
ation. The freer type games are more useful for training 
and investigation of the gross effects of changing battle 
conditions on an operation. 
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SECTION 2 



War gaming in any form may be characterized as a time 
varying process in which decisions must be made in the course 
of a conflict of some type by the opposite players. The vari- 
ables which describe the game (the state variables) change in 
a framework of time called game time, which may be in discrete 
increments (every half hour, day or week, etc.) or may be 
continuously varying. The decisions which are made within 
the framework of the rules of a particular game have an 
effect on the forces of each side which would be instan- 
taneous, if some time delay were not provided for in the 
rules of the game, and if the model is to have any foundation 
in a reality which is determined by time and space, so must 
the model reflect action in time. 

In a game, these time delays for the occurrence of an 
event may be assessed against a side in three general ways. 
First, for actions of a particular type in a particular situ- 
ation a fixed time delay , always the same, may be assessed. 

This time may be multiplied by a degradation factor for weather, 
light, climate, etc., but the base time will remain the same. 
-Second, for actions of a particular type in a particular situ- 
ation, a random time delay may be assessed. This time delay 
will be determined from a known or assumed probability density 
function which ideally approximates the probability density 
function of the real event. Thus we avoid the pitfalls of 



having an always perfectly functioning military machine 
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operating in the game and, while more difficult to handle than 
fixed delays, the extra effort may be justified by the type 
investigation that the game is intended to make. Third, 
arbitrary time delays may be assessed which arise out of 
factors outside the game itself. With no known or assumed 
underlying probability density function, they are products of 
expedience (to speed up closure of forces, for instance). 

Fixed time delays are constant delays in game time as- 
sessed against a particular event. These delays are in 
general an aggregate of specific time delays set down in a 
rule book and held constant throughout the game. 

Fixed time delays have some distinct advantages. In a 
mechanical game they are easy to program and require a mini- 
mum amount of computer core memory in the game iterations. 

In manual war gaming the use of fixed time delays 
greatly simplifies the already complicated assessment process, 
and materially aides in reducing the time required for a play 
of the game. 

Fixed time delays can frequently be obtained from ex- 
isting manuals which give a time for various events. In 
cases where there is no published estimate for the time to 
complete an event the assessor group or the programmer must 
establish these times prior to the play of the game. 

There are obvious disadvantages to the use of fixed 
time delays. In war gaming, fixed time delays allow the 
thinking antagonist an advantage in coordination which he 
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would seldom be afforded in a real life situation. As we 
will subsequently demonstrate, there are instances in which 
the use of fixed time delays can well raise questions as to 
the validity of conclusions arrived at through simulations 
using this type of delay. 

Arbitrary time delays are adjustments in game time to 
maintain the flow of the game. There are obviously situ- 
ations in which rigorous application of fixed time delays in 
game time may produce impractical delays in the real time 
progress of a specific game while yielding no contribution 
to the basic objectives of the game. In these situations 
the control group or programmer may assess arbitrary time 
delays; or no time delay at all, for some event, or events, 
in order to maintain the flow of the simulation at a practi- 
cal rate. The advantage and necessity of this procedure is 
undeniable; however, it must be recognized that such action 
is a reorientation of the game and as such must be considered 
as a subsequent phase of the original game. 

Random time delays are random variables obeying an as- 
sumed probability law used to better approximate real time 
in a simulated environment. The use of random time delays 
will in many cases lend validity to a war game in that it 
denies the undue advantage to the planner of exact coordi- 
nation of rendezvous times, fuel and ammunition requirements 
as functions of time etc. 
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The use of random time delays overcomes many of the 
objections to the use of fixed time delays in simulation, 
however, compared to the fixed or arbitrary assessment of 
time delays the method of random delays is slow and cumber- 
some. 

Since there are few instantaneous reactions possible 
on the part of a military commander, time delays to reflect 
normal reaction times are very important in war gaming simu- 
lations. These time delays are generated in response to 
specific game actions and vary in length from minutes to 
hours depending on the nature of the event being simulated. 

One can easily imagine real life events in which a few 
minutes difference in the time required to perform a spe- 
cific task can produce a difference of magnitude in the out- 
come of an event. In a like manner time delays assessed 
against planned actions in a war game can produce marked 
differences in results. 

As an example, let us consider a simple war game, the 
"Air Battle Analyzer." (Appendix (1)) 

The "Air Battle Analyzer" is a scenario type, free- 
play war game in which each opponent is given a specific 
capability and plan of action. It attempts to provide, and 
display, a readily accessible means for recording and dis- 
playing chronologically the principle movements and operations 
of the different surface and airborne units involved in a 
battle. This display points up the interactions between 
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different units, in particular their ordering in time. 

Let us look at the description of the First Intercept 
Range nomograph, enclosure (1) to Appendix (1). It will be 
noted that a delay of 3 minutes was arbitrarily selected as 
the "delay for decision, " etc, A glance at the Intercept 
Range Versus Detection Range nomograph shows that it is 
scaled for delays from zero to four minutes. 

Using this nomograph, let us look at the intercept 
ranges using delays of l/2 minute and if minutes rather than 
the delay of three minutes assumed in the example. We note 
that an assumed delay of l/2 minute, a detection range of 150 
miles and a target 3 peed of Mach 1.2 yields a detection 
range of 95 miles while an assumed delay of i|. minutes holding 
the detection range and target speed parameters fixed at 150 
miles and Mach 1.2 respectively yields a range of first inter- 
cept of 65 miles . 

Using the Fire Power Analyzer with first intercept ranges 
of 95 miles and 65 miles we obtain the expected number of 
intercepts for 95 miles as 6.9, and for 65 miles, 5»7. 

Entering the target speed scaling nomograph with these ex- 
pected values and assuming a kill probability of .7 for our 
missile while again using a target speed of Mach 1.2, we ulti- 
mately obtain 6.1 as our expected number of kills for a 95 
mile intercept range and 3.9 as the expected number of kills 
for a 65 mile range of first intercept. 

Here one should note that the Analyzer missile system 

is defined as a one launcher, one director system. One can 
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hardly assume that the expected number of kills for a single 
missile system can be multiplied by the number of launchers 
in the entire defense to estimate the grand total of expected 
kills; this would assume independence which does not exist in 
reality. As an illustration, let us assume that we are simu- 
lating 10 missile launchers defending a perimeter against a 
large wave of attacking aircraft. Further, let us assume 
that at least 1/2 of the missiles launched have a target 
which can be engaged by no other missile. This is at least 
the equivalent of f? independent launchers. Then, $ (2.2 
expected kills each) implies a difference of at least 10 
attackers destroyed depending only on the extremes in time 
delays that could be assessed. Were random time delays used, 
discrete times (say fifteen seconds apart) could be randomly 
selected according to a probability law approximating the 
true distribution of delays and allowing a positive proba- 
bility of selection of each discrete time including the end 
points . 

• If the assumed time of l/2 minute or of Ij. minutes was 
known to be correct in every case, there would be no justi- 
fication for the use of random times; however, it is intui- 
tively obvious that such a "correct time all of the time" 
does not exist. In such cases, presumably a random time 
from a distribution approximating the true distribution of 
times for such events would be more realistic, on the average. 
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Thus in the above example, the magnitude of the differ- 
ence in the expected number of kills, and the knowledge that 
the delay for a decision frequently varies with the decision 
maker and the situation, would seem to indicate that the use 
of random time delays was in order. 
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CHAPTER I 



GENERATION OP RANDOM TIMES 

SECTION 1 

Random time delays for a chosen event in a particular 
play of a war game are commonly produced by making a proba- 
bility transformation of a random number from a uniformly 
distributed population to a corresponding number from a 
population which is assumed to follow the true probability 
law governing the event. 

By appropriate scaling, which will be covered in a 
later chapter, this random number is converted to a delay 
time, which is a number within the range of those possible 
for the event. It is this number which is then assessed 
against one of the players as the time necessary for the 
accomplishment of that event. 

Delay time is itself a random variable, that is, a 
function whose possible values lie in some continuous inter- 
val of time. The interval is defined as the range of possi- 
ble values the delay time may take or. 

That the observation of a specific delay time for an 
event is a random phenomenon, i.e., one obeying probabilistic 
rather than deterministic laws, is obvious from any obser- 
vation of the real world. The function which assigns a real 
number to each outcome of the event or experiment is called 
a random variable and is denoted by X, Y, Z, etc., and a 
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particular outcome of the experiment is denoted by x, y, z, 
etc. 



The probability density function, f (x), of a continuous 
random variable X, is a function such that 
f(x) ^ 0 where xeR and R = (a,b) 
f(x) = 0 elsewhere 
and such that 

/f (x)dx = 1 

A probability density function for a continuous random vari- 
able, gives through integration, a description of the rela- 
tive frequencies of the occurrence of particular intervals 
on the real line. 

Associated with these relative frequencies or proba- 
bilities of numbers in a continuum is another function, 

F(x), such that if x is the number associated with some 
outcome of the experiment, then 

F(x) = Probability (X ^ x) = P (X ^ x) 
and, F(x), is related to the probability density function, 
f((x), by 

f(x) = F(x) 



that is 

F(x) 



f (x)dx 



-o<=> X 



and is called the cumulative probability or distribution 
function. (1^) 

If Y is a uniformly distributed random variable on the 
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interval Jo, i] , then 

r (y) = 1 

= 0 elsewhere 

Letting Y = G(x), where G(x) is the cumulative distribution 
function of the random variable, X, which we may identify as 
the true distribution of the random variable, X, the delay 
timej then by letting 

x = G'^y) 

we have a correspondence between the random variables X and 
Y and, through a knowledge of the distribution function of X, 
namely G(x), are able to pass from a knowledge of a particular 
x to a particular y and from y to x. If G(x) is the true 
distribution of the random delay times of an event, we may, 
by taking a random number from a uniformly distributed popu- 
lation, produce a number from the true distribution by the 
transformation above, called the probability integral 
transformation. (If?) 

This generation of the random time delay for the event 
presupposes the knowledge of the true distribution of the 
times for the accomplishment of the event, that is, 
knowledge of G(x). Indeed with such knowledge, the gener- 
ation of random time delays would be relatively easy, albeit, 
laborious. 

In practice, the true distribution of times for the 
accomplishment of the event will generally be unknown and. 
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without this knowledge, we can proceed no further with the 
generation of random time delays and must perforce resort to 
arbitrary or fixed time delays. 

If we would persist in the desire for the introduction 
of randomness in the war game we still, however, have two 
courses of action open to us. 

These are: First, assume some distribution of delay 

times arbitrarily; for instance, the uniform distribution 
might be selected or, for that matter, any other distribution. 

Second, seek some knowledge of the distribution of the 
particular event of interest in the real world. 

The remainder of this paper will investigate the second 
course of action. 

The true distribution of the times for the accomplishment 
of the event can be approximated with some knowledge of the 
times that the event actually takes when occurring in the 
real world. These times (the observations), if they are a 
random sample from the true distribution, are independent. 
Furthermore, if they are dependent in a known way, we may 
still apply many standard statistical techniques. 

Having such data, we are faced with a standard curve 
fitting problem. First, to determine a curve which best 
approximates the true distribution and then to determine 
some quantitative measure of its goodness of fit. This 
measure reflects the "confidence " with which we will state 
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that the curve approximates the true distribution. Classi- 
cal curve fitting consists in estimating the parameters of 
some distribution, assumed to be the true distribution on 
intuitive grounds, by using moments calculated from the 
sample. A further refinement, called the method of moments, 
developed by Karl Pearson in a series of papers (22) from 
1895-1916, sets up a system of objective criteria by which 
a particular curve would be selected to approximate the true 
distribution. The selection of a particular curve depends 
on the values of certain quantities which are functions of 
the first four moments of the sample and the system presents 
a method of calculation of the curve. 

The Pearsonian sjstem is essentially the solution of 
the differential equation 

d£ = y(m - x) 

^ (a+bx+cx 2 ) 

where the values of the parameters, a, b, c and m, determine 
the shape of the curve for a particular treatment. For ex- 
ample, the normal curve is generated if b and c are zero, m 
is the mean and a is the variance. Under this system, the 
sample mean m, and the sample moments about the mean, 
m^, and m^, (i.e., the first three central moments) are 
calculated and are used to produce certain statistics de- 
pendent on the sample size, N, (the number of observations). 
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These statistics are : 



k = ra 

1 

k 2 = 



k 3 



V 



N 



'('N - 1 ) (N 2) a 2 
N 2 

(N - 1)(N - 2) m 3 



Tn 



n 2 Rk+d^ - 3(tl - 1 1 * h ] 



(N - 2)(N - 3) 



where s is the sample standard deviation 



4= Jfc 



These statistics are used to provide estimates 

of and ^ which are Pearsonian measures of skewness 
and kurtosis. These quantities ^ and 4 are used to calculate 
.4’inr and ^ ~t~ 3 ? which are tabled to indicate the 

proper form of the curve to be used. 



b = - in = 



Further, solution of the relations 

c - — sL 



za+j) 



aa+s) 



a=l-3c 



^ 1 # 

will yield the values of the parameters of the differential 
equation, which need only be integrated to give the density 
function. 

Besides the labor of calculating the moments of the sample, 
there are some problems with the system in that except for the 
Poisson, normal and binomial distributions (18) the moments 
of the sample used to estimate the corresponding moments of 
the assumed theoretical distribution do not possess the mini- 
mum sampling variance; that is, they are not most efficient 
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estimates in the statistical sense. 

This difficulty may be overcome by the use of maximum 
likelihood estimators (19), which, however, require that the 
form of the density function (Normal, Beta, Gamma, etc.) be 
known or at least assumed. This is not entirely unreasonable. 

Basically, the principle of maximum likelihood is to 
consider every possible value that the true parameter of the 
distribution, which we have as sumed, might, have and for each 
of these values compute the probability that the particular 
set of observations at hand (the sample) would have occurred 
if that value were the true value of the parameter of our 
assumed distribution. 

Now of all these possible values of the parameter, we 
would choose the one for which the probability of our sample 
occurring would be the greatest. Instead of working with 
the probability density function of the assumed distribution, 
we may work with its logarithm L (purely as a matter of con- 
venience). By elementary calculus if 

dL =0 and d 2 L x _ 

35 * 35 “^° 

where 0 is a parameter of the assumed density function; we 

indeed have a maximum and solution of the equation, 

dL = 0 

35 

for ©, will yield a maximum likelihood estimator for 0, 
called 0. In the case of the maximum likelihood estimator 
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of 0, if it exists, it can be proved that the quantity 
107 ( 0 - 0 ') is asymptotically normal; and that 0 has 
the smallest limiting variance of any estimator. This im- 
plies that, as N becomes very large, $, the maximum likeli- 
hood estimator, becomes most efficient . (20) 

Unfortunately, for this sophistication of the estimators, 
we have sacrificed the crude^ laborious but practical 
Pearsonian System, which has the advantage of requiring no 
prior knowledge of the form of the distribution function for 
its application. But often, there will exist some body of 
past experience which will lead us to choose some particular 
shape or type of curve as being most probable (as, for ex- 
ample, the general concensus that distributions of times 
for simple actions performed by human beings will generally 
be unimodal ) . 

This knowledge or the willingness in a particular situ- 
ation to make an intuitive leap will lead us to assume some 
particular type of distribution to be the best approximation 
to the true one. In this case, the most fruitful method 
will be to estimate the parameters of the assumed distri- 
bution by the method of maximum likelihood, assuming such 
estimates exist. These methods or derivative ones such as 
the Gram-Chalier Series, Edgeworth Series, etc. (18) all 
provide the equations of curves which will fit the data to 
a greater or lesser extent. 
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To recapitulate; having been given a set of data, which 
can reasonably be assumed to be a random sample for a par- 
ticular event, we would apply the Pearsonian procedure to 
determine the shape of the curve to be used as the approxi- 
mation to the true distribution in the case of no prior 
knowledge about the nature of the true distribution. If 
some knowledge of the form of the true distribution is 
available, we would estimate the parameters of the assumed 
distribution by the method of maximum likelihood. The 
principal methods of establishing a quantitative measure as 
to what extent the particular distributions that have been 
derived or assumed to be true by the above approaches actu- 
ally approximate the true distribution are: the Chi-square 

Test and the Kolmogorov Test. 

The Chi-square Test ( ~)C) is a measure of the degree 
to which a series of observed frequencies deviate from as- 
sumed theoretical frequencies over the whole range of values. 
The relative magnitude of this discrepancy is defined by 
the quantity 



where f is the relative frequency of the observations in 
the interval t, f^ is the relative frequency of the assumed 
distribution in the interval t, n, is the number of intervals. 
This value of 'jC' is then compared with the tabulated values 
of the X distribution which give, for various degrees of 
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freedom, the probability with which one would expect such a 
value of X . We may therefore set up a test of the hy- 
pothesis that this sample indeed comes from some assumed 
distribution by calculating the value of X and, on comparison 
with the tabulated values, accept or reject the hypothesis 
depending on whether the observed value of X is greater or 
less than the theoretical value at any probability level we 
desire. Now it will sometimes occur that several assumed 
distributions will fit the data and, in this case, there is 
some justification for choosing the one with the minimum value 

-v/ V 

of K as the curve of best fit. But extreme caution must be 
exercised in the application of the Chi-square Test. The 
reason for this is that no account is taken in the test of 
the distribution of the discrepancies between f Q ^ and f^ at 
particular intervals over the range. It is conceivable, 
because of the nature of the statistic, that, having calcu- 
lated a value of Chi-square by fitting a curve very closely 
over part of the range and relatively far away in the other 
part of the range, this value would be lower than that which 
would be produced by calculating a value of Chi-square for 
a curve which displayed a more constant discrepance over the 
particular intervals. Yet this distribution of discrepancies 
may materially influence our judgment as to the goodness of 
fit. 

For small sample sizes, where the Chi-square Test 

breaks down, we have still available to us the Kolmogorov- 

Smimov statistic, which, in essence, measures the maximum 
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vertical difference between the sample cumulative distri- 
bution function (a step function) and the assumed distri- 
bution. Then comparison of the value of the statistic is 
made with the tabulated value as in the Chi-square Test and 
acceptance or rejection made on the basis of the magnitude 
of the computed statistic. 

The statistic (21) for the two sided case is; 

D n = SU P l p n ( x ) " p ( x )l 
n all x 

where F(x) is the assumed distribution function 

and 

F(x) =0, x ^ 

“ £ X U)* xiX ( j+ i)’ 3 = 1 '— > n ‘ 1 
■ 1 > x - X (n) 

is the sample cumulative distribution function of the ordered 
observations . 

The statistic is itself the least upper bound (sup) of 
the difference between the cumulative distribution functions 
taken over all the values of the argument. It should be 
noted that whereas the assumed distribution F^, is con- 
tinuous and non-decreasing, the sample distribution p n ^ x j* 
will be discrete and so constant over each interval, for ex- 
ample. 
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Di(P) can be considered as having come from a random 
sample on a uniform distribution in the interval £o,l3 (19) 
and, hence is independent of the choice of P, so the 
Kolmogorov-Smirnov Statistic is distribution free. This 
statistic has an advantage over the Chi-square Test in that 
it deals directly with the value of the observations and 
avoids the problems inherent in the summing process of the 
Chi-square Test, which, as has been mentioned, may tend to 
obscure lopsided fits. The choice of one or both of these 
tests of goodness of fit will, of course, depend on the 
particular situation and the degree of approximation desired. 

This brief summary has been intended as the general 
outline of the procedure to be followed, if sufficient data 
is available for approximating the distribution of delay 
times. It should be borne in mind that the importance of the 
knowledge of an approximation to the true distribution 
function must justify the expenditure of time and effort 
necessary to obtain it. 
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SECTION 2 



In order to apply the techniques discussed in Section 1, 
it is necessary that the analyst be provided data. Pre- 
sumably the procedure for obtaining the requisite data would 
be to analyze times measured for particular events under 
actual combat situations. We can immediately foresee diffi- 
culties in this approach. 

First, there is little data recorded in sufficient 
accuracy and detail to provide adequate estimators for such 
parameters as the analyst might deem necessary. 

Second, were the data for some particular event availa- 
ble, one could hardly assume that this data would yield the 
average times associated with similar events unless the 
event had been run, and the times recorded, a large number 
of times. 

Last, the time required and the expense of the neces- 
sary data search is prohibitive. 

One might next logically consider conducting experi- 
ments in real time situations and simply recording the 
desired times in the context of the experiment. Quite a 
bit of work has been done in this area. 

In the course of this thesis, a careful search was 
conducted of existing published data. Contact was es- 
tablished with the U.S.A. Combat Development Experimentation 
Center, Fort Ord, California, and a thorough search of their 
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published and unpublished experiments yielded little of value 
in determining the underlying distribution of the event being 
simulated. This is not surprising. The experiments conducted 
by CDEC, like those conducted by similar organizations, are 
extremely complicated. In each case, the experiment is de- 
signed to consider only some specific portion or element of 
the entire event. The use of scenario type situations and 
arbitrary time delays, as defined in the introduction to 
this thesis, preclude the results of the experiment being 
used to describe the distribution of the entire event. 

In general, the conduct of experiments in real time 
situations to estimate the parameters of an event in an 
actual combat situation is prohibitive for one of the three 
following reasons : 

First, the actual combat situations desired are almost 
impossible to duplicate. 

Second, the near simulation of actual combat situations 
is, not only dangerous to the participants, but prohibitively 
expensive. 

Third, the number and diversity of the events it is 
necessary to describe is beyond effective description. 

In light of the above, we are forced to conclude that 
in general, the procurement of the necessary data, from past 
actions, or by experimentation, is not practical for the 
majority of military events. 



2k 



The following chapter will suggest a third source of 
information by which it is felt that a form of data suf- 
ficient for the description of the times associated with 
military events can be obtained. 
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CHAPTER II 



AN ALTERNATIVE USING MINERAL DATA 

SECTION 1 

In Chapter I we discussed the difficulties entailed in 
attempting to describe military events using recorded or 
experimentally obtained data. 

There is a third possible method of arriving at the 
desired parameters for the underlying time distributions we 
desire to describe. This method entails the use of pro- 
fessional estimates. Professionals, utilizing their 
background and experience in a particular field . should logi- 
cally be able to give a reasonable estimate of the time to 
complete a specific event under certain conditions. 

This method offers advantages over the two mentioned in 
Chapter I in that it is cheap, and with the availability of 
professionals in the field of interest, quite rapid. It 
suffers from two obvious shortcomings : 

First, what estimates are you going to ask the pro- 
fessionals to make? 

Second, how valid will these estimates be? 

There is really no good way to check on the validity of 
subjective estimates short of running a series of the same 
experiment and comparing the mean time of the experiment with 
the subjectively estimated time for the same event. 

The answer to neither of these questions is readily apparent. 
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It would logically seem that in answer to the question, 
"What estimates do you want from a professional? " One might 
reply, "What estimates can he be expected to make with 
reasonable consistency?" Obviously, if we are to request an 
individual to furnish answers based on his experience, we 
must ask questions he has frequently asked and answered for 
himself. 

In this light it is believed that the professional 
military man can interpolate or extrapolate from the domain 
of his experience to furnish the following four estimates 
for many events. 

1. The shortest possible time in which the event can 
be completed, call it 

2. The most likely time required to complete the event, 

say t m . 

3. The most pessimistic time to complete the event, 

A 

tp . 

■ a 

q. The percentage of times, P, the event takes place 
within some given time, centered about t . 

Where t . % , and $ are estimates of the real values 

o' m p 

to, tm, tp and P. 

Before we attempt to answer the second question, let us 
examine subjective estimates in greater detail. Subjective 
estimates are not new. The military commander always has, 
and probably always will be required to use his estimates 
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and those of his staff to effectively plan an operation. 

How does he check his estimates? If the operation comes off 
somewhat as planned, he has done xtfhat we are unable to do. 

He has estimated the time for an event and then run an ex- 
periment which furnishes an immediate check on his estimate 
of this time. Prom each operation he learns. It is this 
equity in experience from which we must draw our estimates 
for events we cannot otherwise practically measure. One 
xrould generally feel that if a time estimate from one 
lmowl edge able individual was good, that the average of the 
estimates of ten of these professionals would be better. 
Ideally it would be hoped that as the number of qualified 
persons estimating the time for some event was increased, 
the average of these estimates would approach the time value 
of the parameter being estimated. 

It is under the above assumption that the use of sub- 
jective estimates is justified. 

How are subjective estimates obtained? The pollsters 
have long realized that the obtaining of subjective estimates 
from a population is a delicate and difficult problem. 

Through experience, they have come to realize that the phrasing 
of the question, the occupation of the individual, and even 
the time of day can affect the answer received to a particular 
question. It may well happen that a wide variation in the 
response to a particular question is the result of a poorly 
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defined or ambiguous question rather than a demonstration of 
the recipients inability to furnish the hoped for response. 

In general, subjective estimates are obtained through 
either an interview or through response to a written question- 
naire. Both methods have their shortcomings. The written 
questionnaire generally suffers from poor response, misin- 
terpretations due to ambiguous x^ording, and the failure to 
distinguish between those persons unable to answer the 
question and those who are able, but unwilling to answer. The 
interview (while probably furnishing more accurate information) 
is more costly, frequently inconsistent due to variations in 
the delivery of the question, and difficult to quantify. 

Since in our particular case point estimates of specific 
times are desired, the use of written questionnaires would 
appear to be the more practical approach. 

In an attempt to overcome some of the objections to the 
written questionnaire, the following rules to be used in 
obtaining subjective time estimates of military events have 
been set down. 

1. Select a group of military professionals capable of 

furnishing time estimates for the requested events. 

It is worth noting, that since we are after ^ ^ , 

O m 

A A 

t and P for particular events we must choose the 
group to fit the events, not the events to fit the 
group. Where possible, a homogeneous group comprised 
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of individuals with various specialties should be 
chosen. This will allow a comparison of the 
stratification within the group, to guard against 
the consistent, but biased, answer. 

2. Assure response to the questionnaire by follow-up. 
Many persons who fail to respond to a questionnaire 
through indifference are capable of producing an 
acceptable answer. 

3. The questions to be answered must objective; fill-in 
the-blank type questions will facilitate the ulti- 
mate reduction and comparison of data. 

4. The information desired must be defined carefully 
and precisely on the first page of each question- 
naire. Only an introduction should be required 
from the administrator of the questionnaire. Every 
effort must be made to keep the questionnaire as 
short as possible. 

5. The questionnaire should be carefully pretested by 
a trial group to check the adequacy of event de- 
scriptions and guard against ambiguities, etc., 
prior to general use. 

6. The data resulting from the experiment must be 
carefully evaluated. 

Following the above listed guide, a questionnaire 
covering ten hypothetical military events has been prepared. 
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(Appendix (2)). This questionnaire is directed to the ex- 
perienced Marine Corps officer. The officers are assumed to 
vary in rank from Captain to Lieutenant Colonel, all with 
service in excess of ten years. 

Let us take the definitions for t , t and t from our 

o m p 

sample questionnaire. This more precise definition will be 
required when techniques are discussed for analyzing the 



data from an experiment. 

t = most optimistic time = The minimum time required to 

accomplish the action if all 
circumstances, personnel re- 
action, weather, terrain, 
etc., combine to work in your 
favor. 



t m s most likely time = The time required to accomplish 

the action under "normal" circum- 
stances. "Normal" to be con- 
strued as those day to day 
conditions you would expect to 
encounter under circumstances 
such that the described action 
could logically be assigned. 



tp= most pessimistic time = 



The time required to accomplish 
the action under extreme adverse 
conditions, and poor response on 
the part of participants. The 
conditions are not so severe as 
to preclude accomplishment of 
the mission, but weather, fatigue, 
human error, etc., all exercise 
strong degrading effects. 



Inspection of the questionnaire will indicate that each 
question is f ollowed by the statement : 

"Percentage (%) of times the event will occur within 
<C hours/minutes of the most likely time." 
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This question is designed to obtain a feel for the a- 
mount of the distribution located within a certain range of 
the most likely time. Answers to this question afford the 
experimenter a method of obtaining a rather crude, but 



nevertheless necessary estimate of the central tendency of 
the distribution. 




2<f 



where the shaded area tinder the arbitrary curve is the per- 
centage estimate of the times an event will occur within <£ 
units of the most likely time. 

It is apparent that the choice of is important to the 
outoome of the experiment in that the ultimate selection of 
an approximating distribution may depend on h. As can be 
seen from the figure, the smaller our choice of cf , the more 
accurate our estimate of h. It is therefore desirable to 
choose the smallest interval (2<f ), for which the experimenter 
can obtain consistent estimates of the percentage of times an 
event will fall in this interval. The question to be answered 
is probably more psychological than statistical. How small 
an interval can be chosen such that reasonably consistent 
estimates can be anticipated? 
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In the questionnaire (Appendix (2)) <S was chosen such 
that 2 £ == R or 25 per cent of the range. 

In most cases the experimenter can himself furnish a 
crude estimate of the range of the event to be considered. 

It would then be a relatively simple matter to prepare sever- 
al questionnaires using the same questions, but varying & as 
a function of the experimenters estimate of the range. In 
this way one could obtain a feel for a subject's ability to 
discriminate between various ranges about the most likely 
time. 



Before we proceed to describe methods for evaluating 
data, it is best that we discuss some of the assumptions we 
are forced to make in conducting our experiment in this 
manner. 

When we decide to use subjective estimates as a means 

of securing t , t , t and P, we are assuming that our 
o m p 

sample values are distributed in some fashion about the true 

values of these parameters t , t , t and p. 

o m p 

V/e further assume that these estimates are so distri- 
buted that the central limit theorem holds. Thus, where N 

A 

is the number of estimates, t Q , the mean of our estimates 
for t Q , is normally distributed with mean t Q and variance 
decreasing to zero as N-* co , irrespective of the underlying 

•A A /\ ^ 

distribution of t^. The same holds for t , t and p. 

o m* p ^ 
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Therefore, under the foregoing assumptions, and for 
fairly large N we feel that we know t Q , t^, t and P quite 
accurately. 

We are fortunate that our method of sampling by question- 
naire is relatively inexpensive. This allows us to choose a 
sample size large enough to insure that our experiment falls 
within the domain of large sampling statistics as defined 
above. This permits us to consider our estimates for t Q , 
etc,, as being approximately normally distributed. 

Let us assume that 100 of the questionnaires in (Appendix 
(2)) have been duly completed and returned. What will we do 
with the data? How will we decide whether these estimates 
are adequate for our needs? What tests should we run? 

Vie want estimates for t , t , t and h. Let us first 

o m p 

compute these for each event. 



100 

? t t •• 

L o IOO c* ''oy.* 



t -i 



A 

where t 



olj 



l- ) J i, *•*, IOO 

J - L 2, • • v lo 

is the most optimistic time 



estimate from the ith questionnaire for the jth event, 

A 

t and t . would be calculated in a like manner, 

Pj mj - 

A A 

As previously defined, h = Percentage , Thus, h; = _fc2_ 
i co 

= 1 Vn. • Where p. * ia the estimate from the ith 

355* Z p °' u 

questionnaire for the percentage of times the jth event would 
fall within t £ of the most likely time. 

Taking say the jth event, let us make a linear plot for 

A A 

t . and Py , to investigate 
Pj 



A 



A 



each of our parameters t Q j ^ t^ 



m J } 



% 







r 



the possibility of outliers or sample stratification. 

When we take our sample of the time estimate for a 
particular event we are assuming that all of our observations 
in the set, comprised of our sample constitute a random 
sample from the same population. For various reasons there 
will occasionally occur in the sample a small number of 
observations which differ considerably from the others. 

These observations may cause us to doubt that all of our 
observations are truly from the population we Intended to 
s ample • 

A 

As an example let us assume we have plotted t for the 

jth event. Let us further assume that this linear plot indi- 

a 

cated that 99 of our 100 estimates of t Q were reasonably 
grouped, but that the remaining, say the kth observation was 
significantly different from the rest. We cannot automatically 
assume that this kth observation is an outlier and thus delete 
it from the sample. 

We can however delete it for cause. If on checking the 
kth questionnaire we noted that it had been completed by a 
twenty-one year old, 2nd Lieutenant, Disbursing officer, 
while the event being estimated was described as, say, the 
time to conduct some facet of an air operation, we would 
probably be safe in asstiming that our kth observation was an 
outlier. In this event we could delete the kth question and 
recompute our estimates based on a reduced sample N, where in 
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this case N = 100 - (samples rejected as outliers) 



In addition to the rejection of outliers for subjective 
cause there are several objective statistical tests designed 
to reject outliers at a certain confidence level. Since we 
are dealing with a rather large sample, the number of outliers 
that would be rejected statistically that could not be re- 
jected for subjective cause is considered not to be sta- 
tistically significant. 

A 

Again consider the linear plot of say t for our jth 



event. We might well observe that our sample estimates were 
in say two or three distinct groups. In this case it might 
pay to stratify our group by specialty say. Infantry officers. 
Aviators, Supply officers etc. If our event to be estimated 
was described as an infantry action we would probably profit 
by taking our estimate of 1; from the sample of officers most 
closely associated with this type of an event. 

Again this reduction in sample size for cause would re- 
quire a recomputation of all of our estimates. 

For the remainder of the discussion on sample evaluation 
let us define : 



o 



N = the sample sizes as reduced by outliers 
and/or stratification thus for example 



It will be noted that throughout this discussion of 



A 



t 



o 




sample evaluation we have not talked of the accuracy of an 
estimate. As stated earlier, there is really no way we can 



judge the accuracy of an estimate. We can however discuss 
the consistency of our estimates. Probably the best measure 
of the consistency of an estimate, say t m , is the measure of 



When would our sample estimates be rejected? The answer 



to this question will depend entirely on the user. Of the 

ten events described in our questionnaire, very likely some 

will produce variations which would make one hesitate to use 

a 

these estimates, say Var t Q >e , where € is some number es- 
tablished by the user prior to conducting the experiment. 

In this case, the only conclusion we can draw is that either 
the question was ambiguously worded, or that the professionals 
answering our questions were unable to give a satisfactory 
response. 

What can we do if our estimates for a particular event 
have such wide variation as to be considered inconsistent or 
unusable. If the event is of relatively long range and it is 
still felt that random time delays should be used, there 
remains little recourse but to re-word the offending question 
and distribute it to another group of professionals. Where 
the event is of relatively short range and the variations in 
the estimates for the required parameters are large, one might 
feel justified in the use of fixed or arbitrary time delays 
in assessing the time to complete this event. 



its variance. 
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SECTION 2 



We have previously seen that in order to generate the 
random times we seek, we must have some approximation to the 
probability distribution underlying these times. Since data 
of the type usually used to derive such approximations is not 
available, we have suggested use of data of a type which is 
available. In this section, we shall examine what effect 
possession of such data would have on our knowledge of certain 
probability distribution laws. 

Table II-l contains some common probability density 
functions. An examination of the table indicates that in 
each case, f(X) is dependent upon some constants listed in 
the column headed "Parameters." In order to specify the 
probability laws best describing real life events, it is 
necessary to know both the form of the probability density 
function, and the value of its parameters. Hence, if we wish 
to generate a random time to simulate the time required to 
conduct a specific operation, we must have knowledge of the 
form of the probability density function best describing the 
distribution of times, and we must be able to determine the 
parameters involved. 

For the present, we will assume that we know the form of 
the probability density function, and that it is one of those 
listed in Table II-l. We shall further assume that we have 
precise information regarding certain aspects of our event, 

30 
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Table 




which will be the conduct in a finite time of some military 
operation. 

Specifically, we shall assume that the following are 
exactly known; 



t Q - The minimum time which could 
be required to complete the 
operation. 

t - The most likely (most frequently 
encountered) time required to 
complete the operation. 



Henceforth, since our 
note it by T, and its 
P (or) = 



t p - The maximum time which could 
p be required to complete the 
operation. 

random variable is time, we shall de- 



density function by f (t). 
< = I f(t)dt 

Jto 



P [t 



Therefore, 

a > t 

o 



where the lower limit of integration is the minimum time. 
Times less than t Q or greater than t have no significance 
for us. If we make a grpahical representation of f(t), it 
would appear approximately as 




where t Q , the minimum time is the zero of our t axis. t m is 
the point at which f(t) has its maximum value. This is denoted 
as the mode of the density function. 

The last item about which we shall assume to have exact 



ko 



knowledge is 



h - the height of the density function 
at its mode. 

With the assumed information, we shall define 

R = t - t as the range 
P o 

R Q = t m - t as the premodal range 

R^ = t - t as the post modal range 

P P m 

R = R + R 
° P 

Recalling that each probability function has associated 
with it certain parameters we shall now proceed to develop 
knowledge of the parameters from the information assumed above. 

The seven density functions from Table II-l will be dealt with 

in turn. 



a. 



clearly 



The Uniform Probability Density Function 

O &/kS <S1^^7 & 

d~ a — Q Q/^C/ ^ 

^ C £) = io — £ fiz 



„o — ^ — £/=> 

= o & / & e e 

The density function has no mode as such, and its height. 



i.e., h, is everywhere constant, and completely determined by 
the range. 



b. The Triangular Probability Density Function 

fM = 



ah 

_ 2 (k-X 



O 4=X< <y 

0 4= owe/ J^6c)~0 &/se*vAes'e 

xst scale t ^ t fs. t 

P 



h (A>- cy) 

Since the range is 0-"X-/z », we must scale t ^ t 44. t 
into that range. If we let 
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X = 
c? =- 



^ I'o , then 

~£o /e a 



3 ince the mode occurs 
at x = a and t = t m 



and & — — /C ^- 

then /'ft t ) — 2jj6^ £ 0 ) 

/z /e 0 



£0 — ^ — tm 



and h is 
also. 



__ 2(/2~€+?o) 

— O elsewhere 






determined by the range for this density function 



c. The Normal Probability Density Function 

~ y2?r d ex/ ° ( x ~~ ^ J 00 

This density function is determined by Af and d , its 
mean and standard deviation. The mean occurs at the mode 
which implies 

^ 



and C t/r>) — ^ — — — 



or 






/ 



V27rh 



The value of f(t) at a distance 4 <4 from the mean is 
approximately 0.0001. Between (- 7 ' - k <d ) and ( ✓y + 4 <5 )> 

to at least four significant places, the area under f(t) is 
unity. 

In this function, t is permitted to range from -00 to 

+ 00 , and we have limited ourselves to consideration of times 

with the finite range t to t . Two situations can exist, 

op 7 
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1. R q and R can both be greater than or equal to 4 <5 
P 



or 



2. one or both are less than 4 <3* 

Prom the preceding discussion of values of f(t) at 4<3 
it may be seen that in the first situation, we can closely 



approximate f(t) by Q y/ & ‘fr'A 2 



using _ 



-VzttA 



and 



— z6 



•/*? 



9 



The second situation definitely implies that although 
the time required to complete the operation may be normally 
distributed, the density function must be cut off at one or 
both ends (truncation) . This leads to the rejection of use 
of the normal density function under these circumstances, 
and to the truncated normal density function which is 
presently to be discussed. 

d. The Gamma Probability Density Function 

— ^ elsewhere 

If we set d f ( 'X ) = 0 and solve for 25 , we find 

ctx 

that the mode occurs at 

% — <X/3 

Since f ( % ) has its origin at zero, we must rescale our 
time so that t 0 corresponds to 0. To accomplish this, let 



^ “ if “ ^6 
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. 













The mode then occurs at 



x = t - / = /e 0 = or <3 



a#d 



- 






to S if < 



CO 



= <9 



elsewhere 



Using the above, and evaluating the function at t m , 




Making use of Stirlings Approximation 
77/ -O (Z7f)^ 77 * * Q- 77 J 

/ / ^ 






(%)*?* e-% fi 






which reduces to 



/j ^ — or ' y ~ / - 

Y2#Y&g z7r/?o4* 

and (X = ^ ^ ^ / ^<> ^ 

Thus if our maximal time had no bound, h and R 0 are 
sufficient to determine approximations to (X and /3 , 

It should be noted that the approximations depend on t Q , 
t and h, and are not affected by t p . 

Under our conditions however, we do have a finite maxi- 
mum time. This implies that our density function is truncated 
at that point, and in order that 



/ = 




we must adopt 



some alternative. There are three readily apparent from 
which to choose 

(a). l-/e may consider that there exists at t a discrete 

probability mass point where the probability concen- 
trated at tp is defined by 



then 



J <*/ /3 ° ( + / ^ £ ° c/t 

= Oi 7jjpi+i t <t<tp 



l 



0(1 /3* H (t-t.) e~ & 



fa 

- o 

(b ) . we may define 

t} = Aa/rf*" ^ *"*/*<? ~ ^ *„£t < i P 



elsewhere 

£z£* 

A ex* 

— o elsewhere 

t p 



/ tp 

— /u -4 \°c~ - £. 

and is regarded as a constant 



(c). f(t) 



/ 



o<'e 



C*H 



tt<) 



v- 



t- to 



& 



may 



be asymptotically zero at tp, and the probability 
J ^ e/t! = O y ±n which 



case, we regard t p as being so large that its mathemical 
effect is that of co . 






, -»V- 



•- 



















































Under alternative (a), we have in effect established 



that the probability at t may be disproportionately large, 

P 



but have not invalidated our determination of an approxi- 
mation to the parameters o( and /3 , 

Under alternative (b), the probability of a time greater 



but it has presented us with a third parameter in f(t), The 
effect of this is to invalidate our approximation to & and 
/3 . The mode will still be located at £-£ a — && , but 
with only the information at our disposal, we can no longer 
obtain a closed expression for cX , /<3 and A . Numerical 
methods with the assistance of the Tables of the Incomplete 
Gamma Function (16) must be applied to evaluate the 
parameters , 

Under alternative (c), we recognize that there exists 
some finite 



the measure of the probability that a time will exceed the 
maximum, but regard ^ as being so small that it may be neg- 
lected, This implies that 



than tp has been redistributed over the interval 





where ^ is 




O 
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Since 0{ > - / and /3 > O , the constant tern 



can only have the effect of reducing f , and hence, if 

/Z°‘e-'^=o J =>/Yzt ,)± 0 

This will occur if 

» /e* or 
^ » <* /€ 

but Ctv3 = /fo 

so that if 

/? 



.4,/e 



» /e* , we may disregard T? , 



e. The Beta Probability Density Function 

(c ) ■f('x)= o£X< i 

Of 0 /O i 



= o 

The mode occurs at 



elsewhere 



■X = * & 

Since the range is 0 - 1, we must scale our times as 
follows 

Let ;* = t. Z*1 - t - go 



- to 

* = ^ 



/e 



The mode occurs at 4 , or % - 

/. = /e 0 



_&o 



=* % = % - * 

( u ) a = JL/Ze 

fS = ^ where k is some constant 

to be determined. 

Evaluating ^' ) at R& 

"'fit. * 

A = -L t 'll 

” * " a*)/ w('-«) 

where the factor -J- is required by the scaling . 

fa 

Simplifying * - 

. /&r 

a&xufrM { *' ' 

Again making use of Stirlings Approximation, 

a*H)vznU *.)***■%-** 

+**<?** 






which reduces to 

fa/? ^ 



faH 



/e 



/ Z7T V/ifa* fas 



or 



(tot) - = A -iAtt fa fa fa P 

a relationship which can be solved to provide a value of 
h. Since from (it ) 0( = <Jz fa* 

ft — ^fa/ 9 ^ the solution to 

C ‘‘‘ ) completely determines Q( t/3 in terms of fa* € fa/* 
it should be noted, that the equation {tit ) in solution 



places a bound upon h. 



fao fa/o 



as lower values 
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of h would give imagirary roots in solving for k . 

This implies that although /?<> and do not determine 
h for the Beta distribution, they do provide a lower 
bound, 

f. The Log Normal Probability Density Function 



ySTdx exp ( } x> 



= o 



elsewhere 



The mode of the density function found by setting 
~f(~x) = o and solving for x is 



x =■ e 



— CS ‘ 



Postponing discussion of the upper limit of the range 
we scale 

so that t - to ~ X ' 

,2 



At 



^ Z/w to j /to = Q 



'V-d 



or 



/r/e ° ) = 



or V27rA /£> = 



ex/o 



f-*7 



z rA = t 3 - e d 



upon squaring 



a value of n 



3y means of a simple graph of & 
can be selected such that 77 & _77 = ZTfr? /?o . 7/77 

will then be <3 and 77 will equal . 

As was the case with the Gamma density function, we have 
so far been considering a density function whose variate is 
without upper limit, as applied to a time variate with a 
finite upper bound. Also as was the previous case, our de- 
termination of the parameters was independent of this upper 
bound. 



In general, we may choose from three alternatives to 
resolve this problem. 

a. Regard t^ as a point possessing the following 
probability mass 

(/) 1 Vzvd(d-te) eX/D '*) } ** 



which leads to 



for 

~ (0 above for t = t 



elsewhere 



As long as t > t m , the method of evaluating ^ and d 
remains valid . 
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/ 



b. Regard 

fU)= 



AV2r6(t 

for 



— exp 



t 0 * 1 f= fp 



— O elsewhere 

where W = I V&rd(t-&J 

Cq ' 

This has the effect of introducing an additional 

parameter A to consideration, and imposes the 

requirement for at least one item of information other 

than t , t , t and h for a unique determination of 
o m p 

'V/ <5 and A. 



m) -J] 



may be asymptotically zero at t 



that 



*P 



which would imply 

O 



t as being so large that the probability of exceeding 

Jr 

it is zero. This requires that 



/ ^ 

~Q J 2, ^ 



(%* y? 1 which reduces to 

/£? >^> dSZ 

5i 



g, The Truncated Normal Distribution 
fCx) = 



A =f SS h *** 



Q<X < h 

•f (x) = O elsewhere 

The mode occurs at , and our assumed information 



leads to 



vd<SA ex* } 

A ~J expf’&d'tm) ]c/t 



for 



4 Hit* 



■f C t) = o 



elsewhere 



It should be noted that A plays the part of a constant, 
but it is another parameter which must be evaluated. 

Now, 

fftm) = ^ = ^j^7 isa 

constraint on the two pararaeters A and , but with no more 
than this, we could find any number of A’s and <5 s and thus 
any number of possible distribution functions. 

What is required is one further piece of information 
restricting a point on the curve traced by f(t). The value 
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of f(t), the height of the density function at either end 



point (t or t ) would be sufficient, 
o p 

If such information were in our possession, it could 



oe used as follows : 



/ 



\ ' A/r „ d 






and from 
which leads to 



i) 






_ 

Q 2d 2 ’ 



ii) <d 






could uniquely determine 

If however, all that were known was that f (t) was 

asymptotically aero at t , then the above would provide no 

P 3h 

assistance in determining ^ . Under the assumption thougli 
that f(t) = 0 at t , and that (f(t) 0 at t 0 (i.e., R a <, 

dp ), then .5 ^ “• • This follows since t ^ t Q and 

tre have not truncated more than l/2 the area under the basic 
normal distribution . Further, the basic normal distribution 

£(*) = e/ /° (~ (t ■„! 



will 



approach zero sooner than f(t) since 

m) = £zi 'll 

A 



ana 



A <1 



implies f(t) f (t) for any value of t 



We also know from properties of the normal distribution 



that 



fn") 

f(t) 



•= . 0001 = o 



so that 
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or k 0004 T i*'/ < ^ 
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We can proceed to approximate a value for S by 

assuming that f(t) = 0 at t implies 

P 

. OOO 1 < m\ t < , OOI 
Under this assumption, using tables of the normal 

distribution it can be inferred that 

< (t r -t m )< 4<s 

■< /Z P < <4" <4 y and that 

R can be used to obtain an approximation to 
P 



. 0002 



54 



SECTION 3 

Section 1 outlines the general statistical procedures 
to be followed in attempting the normal pattern to determine 
the probability laws underlying a real-life event. Under 
many circumstances, these procedures cannot be followed. 

It may be that the cost of carrying out the necessary repe- 
titions of the event is prohibitive, or it may be that the 
event is an infrequent occurrence of nature, happening too 
rarely to obtain a sufficient sample. At any rate, the data 
necessary to follow standard procedures is not available. 

In the large majority of cases, the times required for the 
completion of various military operations fall into this 

vV; 

category. 

What we shall next attempt is to examine some procedures 
which might be indicated in such cases, but where certain 
specific minimal data is available. 

Given that we have 

an estimate of t Q , the minimum time 
required to conduct the operation, 

/ t_ an estimate of t„., the most likely 
time required, 

t an estimate of t . the maximum 
P P 

time required, and 

/\ 

h an estimate of h, the height of the 
underlying density function at t m . 
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Wh$t knowledge can we deduce of the underlying proba- 
bility law? 

In Section 2 of this Chapter, it was shown how t Q , t m , 
tp and h could be used to determine the parameters of a set 
of Probability Density Functions. It should be recognized 
that at that point, the discussion was conducted from the 
position that there was no search for the form of the density 
function, and that t 0 , t , t and h were regarded as known 
constants of the problem. Here we are dealing with estimates 
of the values of those constants, and have no real knowledge 
of the form of the density function. 

If for the moment, we assume that we have accomplished 
one half our task, and have determined the form of the density 
function which we are seeking, we then have only to estimate 
Its parameters to complete our approximation to the proba- 
bility law underlying the event. Since in Section 2, we 
related these parameters to t Q , t , t and h for seven 
probability density functions, if the form with which we are 

/\ /N /X yX 

dealing i3 a member of that set, we can use t , t , t and h 

o m p 

to obtain estimators for the parameters. If the form of the 
function is not from that set, it will be assumed that ma- 
nipulations of the type performed in Section 2 will produce 
some approximation to its parameters. 

Having thus rather summarily disposed of the problem 
of. parameter estimation, we turn to the much more difficult 
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problem of determining the form of the probability density 
function. The basic information we possess is that we have 
estimators for the following points. 




and are attempting to define a unique curve which has its 

origin at passes through the point a, which is at a 

a a a 

height h above t , and terminates at t_.. The curve is- 

HI P 

further restricted in that its highest point must be at a, 
and, that it represent a probability density function f (t). , 
The latter implies 

f(t) — O and 

C tp 

/ = I 

i 

Under certain circumstances, almost any of the seven 
probability density functions discussed in Section 2 will 
generate such a curve. It must be stressed that these seven 
were selected for discussion only because they are among the 
mo3t frequently encountered, and easy to work with. It would 
be more likely that there are an infinity of such curves, 
each with an associated f(t). With no further Information 
available, we cannot, in general hope to cope with any* such 
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function dependent upon three or more parameters. Thus, we 
are forced to restrict our attention to examination of those 
functions for which we possess enough information to enable 
us to proceed. 

At this point, it should be admitted that we may have 
already excluded the true probability density function^if it 
exists. This should not however be a matter of major concern. 
Even in the more conventional statistical problem (Chapter I), 
no claim is made that the ’'true'’ probability density function 
results. What is a result, is a density function which is 
a member of a set of known, in general tabulated, and con- / 
venient density functions, and which "best fits" the data. 

Even that much depends on far more data than is available 
under our present assumptions. 

Since the seven functions in our set are mathematically 
well defined, and in general, well tabulated, we shall seek 
a "best fit" from that set. We shall first proceed to ex- 
amine each of the various functions in turn to see what 
effect the data can have upon their applicability. 

The uniform probability density function is dependent 
only upon its range t p - t Q = R, and therefore having t Q - 

t - R as an estimate of the range, the function would be 
P 

determined. However, the value of h - 1 is also determined 

R 

by R, This would infer that we should reject the possibility 
that the uniform is our function if h is such that it differs 
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considerably from 1 . 

A very much similar situation exists for the triangular 
probability density function where h = 2. If h is consider- 
ably different from 2, we would reject the triangular from 
consideration. 

As developed in Section 2, the Beta probability density 

2 , I 

function also has a restriction in that h 5 ^ 

/x!o /-£/=> 

A 2, 

Therefore, we would reject it from consideration if h were 
less than ^ -4^ . 

Consideration of the normal density function involves 
several factors. To begin with, it is a strictly symmetrical 
distribution. Unless it is to be truncated, a case presently 
to be discussed, the value of h fixes its standard deviation 
O . If this is the case, and either /^ 0 or /^o are less 
than k <6 > then there are contradictions to our assumptions 
in Section 2, Chapter II which rule out its use. 

In Section 2, it was seen that unless certain assumptions 
are made, our assumed information is insufficient to determine 
the parameters of the Truncated normal probability density 
function. Hence, if we consider this density function, it 
must be with the assumption that the truncation is one-sided 
Even after making this assumption, we are left with only an 
approximation to d> where 

<$•>5 6 < fctp <46 or 
2 ~ 

° v5. S > an( * u* 1 ! 633 
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where 



^ / • / 

cS is such that n = ; — ? ~ 

V27f A 6 

* = j yS^S eXP {'^ (t-t™f} , then 

J £ 

Co 

our assvimptions will lead to contradictions with our data. 

The differences from the above described equalities and 
inequalities upon which to base rejection is subject to the 
users latitude. In all cases, lire are now discussing estimators 
rather than the true nature of R 0 , Rp, and h. The magnitude 
of a difference sufficient to cause rejection can be es- 
tablished only upon consideration of the accuracy of the 
estimates and the desired bounds for approximation. 

For the remaining two density functions, given our as- 
sumed estimators, a probability density function can be 
cranked out which will fit the three desired points, but it 
will have a discrete mass point at t . Should the height of 
this mass point exceed h, then we fail to meet the criteria. 

Any further attempts to determine which provides a "best 
fit" run head-on into the discouraging fact that we possess 
no other points with which to measure "fit". At this point, 
our excursions have reached a stage beyond which future 
progress can result only with more data. Under our hy- 
pothesized conditions, we cannot extract more data. 

Since further theoretical progress seems barred, we 
shall now examine the origins of the problem to determine 
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whether the requirements of the problem and the bounds upon 
an acceptable solution can provide instruction as to future 
routes to pursue. 



Basically, we are searching for a method of producing 
random time delays to assess against participants in a war 
game. These time delays represent the time which might be 
required for the accomplishment of the real life event which 
is being simulated. The only tools considered available to 

A A A 

the assessor are a table of t^, t , t for the various 

o’ m p 

allowed events, and a uniform random number generator. This 
latter is a device, which upon activation, will generate a 
random number drawn from a population with a uniform distri- 
bution on the interval [o - 99~\ . Specifically, we are 
attempting to provide for the assessor another tool — one 
which will convert such a random number to a random time T 
where t Q - T - t and where T has a distribution roughly 
corresponding to the distribution of occurrence times for 
the real-life event . 

There are three requirements for any method to be gener- 
ally acceptable. It must be simple, rapid, and a reasonable 
approximation to reality. 

War games, computer, manual or composite are extremely 
complicated models intended to simulate a most complex real- 
life situation. Any complication introduced into a component 
of play of the game must have its effect measured in the over- 
all game atmosphere. The components do not exist in vacuo, 
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they exi3t to fulfill specific requirements for simulation. 
Unless relatively simple, the complication they introduce can 
over-balance the benefit of accuracy of simulation. ■ 

Present day sophisticated war games are lengthy affairs. 

A manual game may take as much as six months for completion. 
Computer simulations are accomplished in minutes or hours, 
but behind the play of the game is a programming effort that 
is measured in man-years. Any refinement of methods or intro- 
duction of new components can ill afford to further extend 
time requirements. 

Most present war games utilize the method of fixed time 
delays, i.e., each event has associated with it a fixed time 
to be assessed for its accomplishment. Those that differ 
from this method generally apply uniformly distributed random 
time delays, with the implication that all times x^ithin the 
range for an event are equally likely. On no more than intu- 
ition, neither seems a very close approximation of reality. 

Any refinement of method should have some closer agreement to 
the way in which real times vary. 

Having thus reviewed what it is we are trying to do, 
and the bounds on our method for doing it, we proceed to ex- 
amine the effect of these bounds upon selection of a distri- 
bution function. 

The uniform distribution is the simplest of our seven. 

The conversion of a uniformly generated random number 
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requires only a scaling to produce a random time 
This ease of conversion makes any use of a uniform distri- 
bution quite rapid. It has the important disadvantage that 
in many circumstances it is not even a reasonable approximation 
to reality. One hardly expects that in very many events the 
minimum or maximum time will be observed with equal likelihood 
to any other time. 

The triangular distribution is the next simplest. It 
possesses several features expediting work with it. The 
density function is completely determined by the extremes 
of the range and the position of the mode. For a given tri- 

-A -\ 

angular distribution of T J t 0 :5 7~ ^ (determined by 

the position of the mode), there is a corresponding distri- 
bution of 0 - X. £ 1 * By means of a table or graph 

of the applicable £ O — 7^ triangular cumulative probability 
distribution function, a uniform random number Q 0-99~] 
can be converted to a triangular random number O-']. 

This can then be scaled to the interval [ £ - £] to pro- 
vide a random time. This should give a closer approximation 
to reality than the uniform distribution. At least, times 
at the extremes of the range would have the least frequency 
of occurrence, and according to the position of t , any 
degree of skewness may be introduced. If we assume that the 
mode will occur at some discrete point [ 0 - .90] we need 
only concern ourselves with 100 triangular distributions to 
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have a set of graphs or tables which can be used to handle 
any case. In fact, onlj’’ 50 graphs will suffice since for 
distributions skewed to the right, we can use a corresponding 
distribution skewed to the left, taking 1 - the triangular 

-A 2 

random number, and scaling it to ^ 7” — • 

The Beta distribution involves a greater number of 
calculations than the previous two distributions. Its 
parameters must be calculated from the data, and for a fixed 

/S. /\ 

set of values of t Q , t m and t , there are a variety of curves 



which must be restricted to one unique curve by use of h. 

That is to say that for a beta density function on o-'] 
and a given mode, there are many curves possible depending 
on the parameters. The distribution, once the parameters 
have been determined is conveniently handled with the as- 
sistance of Tables of the Incomplete Beta Function, (17). 
Again, by use of a table or graph of the cumulative proba- 
bility distribution, a uniform random number [° - 991 ] 
may be converted to a beta number £o- 1 ~^ which can be 
scaled to ~ t ^TJ . Quite different from the preceding 

case is the number of tables or graphs required to handle 
the possible distributions. Since there is no unique distri- 
bution associated with each modal position, there are many 
graphs required for each such position. Intuitively, the 
Beta distribution may more closely approximate reality than 
the two previous distributions. Its height at the mode can 



6k 



vary independently of R Q , R and R, at least as long as 



h* > 



. Depending on its parameters ck and /3 



/Co /?/> 

the frequencies associated with the extremes of the range, 
particularly the maximum end will be less than with a tri- 
angular distribution over the same range. 

Although not simple, the normal distribution is probably 
the best known of our seven. There exists a greater body of 
literature about it, and more readily available tables than 
for any other. Unfortunately, we have no assurance that we 
can use it in other than truncated form, which to some extent 
degrades the utility of existent literature and tables. For 
any event which has a skewed distribution, non-truncated use 
of it would represent a distortion of reality rather than an 
approximation to reality. For that reason, any further dis- 
cussion of it will be continued in the Truncated Normal where 
the normal may be considered a special case. 

The gamma and the lognomal distributions both have 
features which tend to complicate their use for our purposes. 
Both density functions are unbounded at their upper end. As 
was discussed in Section 2, there are three courses available 
in terminating the functions at t , truncation, consideration 
of a discrete probability mass point at t , and treatment of 

ir 

the functions as asymptotically zero at t . If truncation is 

.t' 

adopted, we can no longer make convenient approximation to 

the parameters. If we adopt the idea of a discrete mass point 

at t. , in generating random times, we will produce t more 
P P 
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frequently than some times less than t . The third altema- 

Jr 

tive will, in general, produce the same effect. Since the 
idea of truncation seriously impairs our ability to produce 
the parameters, its use must be rejected. Under either of 
the two others we can proceed as follows. Prepare graphs 
or tables of the cumulative probability distribution function 
on a scale [»->] . Use these to convert our random number 
to a random number on the appropriate £o-f 3 l°6 norma l or 
gamma distribution. Scale this number 



to [?„ - to 



produce T. For both distributions, the parameters depend 
only upon and h. If we consider that the mode will occur 



at some discrete point in fo-*/3 we must still prepare 

R L J 



many graphs for each such point to consider different possi- 
ble values of h. This, of course magnifies considerably the 
effort required for the use of either distribution. Both 
distributions have the disadvantage that they cannot readily 
be applied to any event which would seem to have a symmetri- 
cal distribution. 

The Truncated normal distribution has the advantage 
that its use is facilitated by the extensively tabled normal 
distribution. If we consider for the present only truncations 
on the lower side, with the truncation at 0 and the height of 
the upper tail as approximately 0 at 1, then for the interval 
[ 0 -. 99 ] , only 50 tables or graphs are required to provide 
cumulative probability density functions for the 50 discrete 
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points ~ at which the mode could occur. This is due 

to the fact that at each position, we have but one curve to 
consider. Using these graphs, uniform random numbers can be 
converted to the appropriate truncated normal distribution, 
and scaled for use as T. For truncations at 1, with the 
lower tail height = 0 at 0, the same tables or graphs can be 
used, but the random number to be scaled would be 1 - the 
random number generated by the conversion. 

Grouped in order of simplicity and rapidity in use, our 
seven functions would appear as 

Uniform Most simple and rapid 

Triangular J 

Truncated Normal Intermediate 

(the normal a special case) 

Lognormal 

Gamma >- Least simple and rapid 

Beta J 

Comparison as to ability to approximate reality is not 
so easily done. In the first place what is reality? For 
all but a limited number of the simplest type events, no 
one knows. It is most likely that no one distribution best 
represents all classes of events. One form of distribution 
may best represent certain classes, but other classes would 
require other distributions. Experimentation has not been 
conducted nor is it likely to be conducted so as to compre- 
hensively determine the distributions associated with the 
classes of events under discussion. We can however be sure 
that under certain conditions several distributions may all 
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bo very close approximations . Even in the standard sta- 
tistical problem with adequate sampling this frequently 
occurs. Using "goodness of fit" tests, txro or more, distri- 
butions are found to "fit" sample data. 

Although we have no factual data as to reality, certain 
generalizations and intuitive conjectures can be made from 
observation and experience. In dealing ivith tines required 
to complete a military operation, it seems plausible that 
many times (if not the majority of the time) one is working 
with an event requiring a skewed distribution. That is, the 
minimum tine is closer to the most likely than is the maxi- 
mum. There aro only so many ways in which to expedite the 
operation, but it frequently seems that there are a million 
ways to degrade it. Also intuitive is a feeling that the 
probabilities of times near either extreme of the range 
should be very loiv. That is, at either extreme, f(t) should 
approach, zero asymptotically. 

Excluding the uniform and the triangular, all of our 
density functions if generated by the same t , t , t and h 

r 

tend to plot close together from the mode to the upper end 
point. Erom the lower end point to the mode, the truncated 
normal x^ill, in general, plot higher than the remainder. All 
of them also conform to our previous conjectures to reality, 
with the exception of the truncated normal at the truncation 
end. 
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From the foregoing, we are reduced to selecting a distri- 
bution function on the basis of simplicity, speed in use, 
and general applicability. Those that are simplest and fastest 
are the least close to what we assume as reality. Those most 
approximate to our notions are also the more complex and 
laborious to use. 

The next chapter of this thesis presents an algorithm 
for use in generating the random times desired. In this 
algorithm, we have elected to make use of the Truncated 
Normal Distribution as the most satisfactory compromise of 
speed and simplicity versus approximation to reality. 
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CHAPTER III 



AN ALGORITHM TO GENERATE RANDOM TIMES 

Section 1 

This chapter describes an algorithm for manual use to 
generate random times. The input to the algorithm consists 
of a uniformly distributed random number, and the user»s 

■A A ^ 

estimates of t^, t and t . The output consists of a random 

O HI p 

time T such that 

A A 

a, t_ ^ T t and 

o p 



b* T is derived from a truncated normal distribution, 
where truncation occurs on only one side and 



3.5 = 



Max 






Lv Bp] 



/\ 

with the mode at t . 






If R o = Rp the distribution is normal, = t 



m 

3.5 <i = R 0 = Hp. 

In accomplishing the above, the algorithm performs the 
following functional steps : 

1, Selects an appropriate truncated normal distribution 



scaled to the interval 



[»->] 



2, Converts the uniform random number to a truncated 
normal random number in the selected distribution, 

3. Scales the truncated random number to Qt Q , t "J to 
produce T. 

The necessary elements to accomplish this are : 
a, A Uniform Random Number Generator 



This is a device which upon activation produces a 
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random number drawn from a population which is 
uniformly distributed from 0 to 39. It may consist 
of no more than use of two columns from a table of 



device to produce such numbers. 

b. A Table of Event Estimates 

This is a user furnished table, listing alphabeti- 
cally, by class, etc., the events occurring in 

A A A 

simulation with their associated t , t , and t . 

o m p 

c. Cumulative Probability Density Functions 

These may be either graphs or tables, one for each 
possible case of the truncated normal distribution 
to be considered. Each graph represents a differ- 
ent possible position of the mode for discrete 



are used, there will be little significant change 
from graph to graph. In Section 4* it will be 
shown that the desired degree of approximation may 
be obtained by considering some number less than the 
maximum. 

d. The Selector Scaler 

This is a mechanical means for selecting the 
appropriate cumulative probability distribution 
function graph, and scaling a random number from 



random numbers, or it may be a mechanical/electrical 



points in 




If the maximum (£l) 



that distribution to the interval 
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Construction of a selector-scaler should be within 
the production capability of any reasonably large 
training aids facility. If no selector-scaler is 
used, the necessary instructions for manually 
carrying out its functions will presently be 
explained. The selector-scaler consists of the 
f ollowing : 

(1) The 0-1 Bar - a horizontal bar with di- 
visions from 0-99. 

(2) The Time Bar - a bar mounted by bracket, 
at a convenient angle to the 0-1 Bar, 
with 150 divisions, and such that it 
may slide through the bracket to place 
any desired dividing point coincident 
with the zero of the 0-1 Bar. 

(3) The Converter Grid - a ruled sheet of 
acetate or plexiglass with equidistant 
parallel lines whose spacing is equal to 
that of the divisions of the 0-1 Bar. 

The right most line is termed the Base 
Line and labelled B. Lines on the left 
half are black, and those on the right are 
red. Each line has a numerical designator. 

A sketch of the components appears as Figure III-l. 

The functioning of the algorithm is as follows: 
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When a random time is desired for an event simulation, 

the event is located in the Table of Event Estimates. The 

Time Bar is moved until t on the Time Bar is coincident 

o 

with the zero of the 0-1 Bar. The Converter Grid is placed 
over the bars so that the base line connects the 1 of the 
0-1 Bar with the value of tp on the Time Bar. The grid 
line on the converter passing closest to t is noted. Its 
number represents the number of the Cumulative Probability 
Distribution Function Graph to be used. A uniform random 
number is generated. Entry is made with this number on the 
vertical axis of the indicated graph. At the point where 
a horizontal from this entry point intersects the graph, 
the value on the horizontal scale is read. This value 
represents the truncated normal random number At 

this number on the 0-1 Bar some line on the Converter Grid 
is superimposed (the base line still connects 1 and t ). 

Jr 

The value that this line intersects on the Time Bar represents 
T. 
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SECTION 2 



Instruction s for Use of the Algorithm with 
the Selector Scaler 

When an event occurs in simulation for which a random 
time delay is desired, 

1, Locate the event in the Table of Event Estimates. 

Note t 0 , t m and t . 

P 

2» Slide the Time Bar until t is coincident with zero 

o 

on the 0-1 Bar. 

3. Overlay the Converter Grid so that line B connects 
the end point (1) of the 0-1 Bar with the value 

A 

of t on the Time Bar. 

P 

ij.. Locate t m on the Time Bar, and note the color and 
number designator of the Converter Grid line 
crossing it. 

5. Locate the C.P.D.P. graph corresponding to this 
designator. 

6. Generate a uniform random number. 

7. With this number enter the selected graph on the 
vertical axis. Where a horizontal from the entry 
point intersects the graph, read the value of the 
horizontal axis. If the Converter Grid line crossing 

A 

t was red, subtract this value from 1. If the line 
liras black, use the value as read. 
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8, Locate the number resulting from Step 7 on the 

0-1 Bar. Follow the grid line crossing this point 
back down to the Time Bar, and read the nearest value. 
This is the desired random time. 

Example: Suppose a random time delay is desired for an 

engineer platoon to install a six row mine belt 
across a valley 1500 meters wide. The Table of 
Event Estimates yields 

t Q = 2 hours 

t m = 8 hours 

t = 22 hours 

The number 2 on the Time Bar is placed coincident 
with 0 of the 0-1 Bar. The Converter Grid is 
placed so that line B connects the 1 of the 0-1 
Bar and 22 of the Time Bar. The black line 
crossing t has the designator .30. Figure III-2 
is C.P.D. F. graph #30 . A uniform random number is 
generated, say 57 - At .57 on the vertical axis of 
Figure III-2, a horizontal line intersects the graph 
where the nearest horizontal scale value is .35* 

At 35 on the 0-1 Bar, the grid line is traced back 
to yield 9 on the Time Bar. Since we are using 
hours, 9 hours represents the random time delay to 
be assessed. 
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SECTION 3 



Instruction for Use of the Algorithm 
Without the Selector Scaler 
When an event occurs in simulation for which a random 
time delay is desired, 

1. Locate the event in the Table of Event Estimates. 

2. Compute ^m - ^o = n, - ^m = n^ - 



3. Locate C.P.D.P. graph # sn/n [ /7 s /7 'J‘ 

Generate a uniform random number. 

5* Enter the vertical axis of the selected graph with 
the generated random number. Where a horizontal 
from the entry point intersects the graph, read 
the value of the horizontal axis. If the graph 
# was n^-, subtract this value from 1. If it was 
n, use the value as read. 

6, Calling the number resulting from Step 5> t, 
compute 

T = ( t - t ) t + t_ to determine a random 
p o ° 

time delay. 

Using the example of Section 2, 




Again using graph # 30, Figure III-2 with a uniform 
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random number assumed to be 57 the value of t 
is .35 

T = .3 5(22 - 2) + 2 = 9 hours. 
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SECTION 4 

Conatruction of Cumulative Probability 
IMLatribution Function Graphs 
This algorithm makes use of graphs of the cumulative 
probability density function (C.P,D.F. ) of various truncated 
normal distributions to convert a uniform random number 
(0-99) to a random number from a desired truncated normal 
distribution* For generality and convenience, the truncated 
normal density is considered to be truncated at zero, and to 
have its upper 3*5 <6 point at one. We are actually concerned 
with a random time T x^hich has a similar density function 

A /s 

truncated at t Q and its upper 3*5 <3 point at tp. The mode 
(and so the mean of the underlying normal density) is at 



m* 



This case can be reduced to the 



[o - l] 



interval by 



scaling. Such scaling reduces to a minimum the number of 

A A ^ 

graphs required, since t Q , t ^ and t p could conceivably take 
on any real positive values. In scaling, any point t, where 



A 



t Q _ t — corresponds to a point 



% 



o — / , by 



letting 



z? - 

& -£ 



zf 1 -Jo 
/€* 



Now the mode occurs at t which could occur at any value in 
[£> , ?/»] . It will suffice for us to consider only such 
cases where 

(fm-Z.) £ 1) or 



<30 



i.e., where truncation occurs at t • Gases where t is the 

o o 

. ✓N. 

lower 3*5 o point, and t is the point of truncation can be 

P 

handled using the same set of graphs, but with slightly 
different application. 

If we desire a set of graphs to cover all possible 
cases, we then require a graph for each possible position of 
the mode. In the interval Qo - 1 J , this would imply a 
graph for each possible value of 





where R 



R 



TT IT 0 

Our uniform random numbers can only take on values in 
(0-99)# and the truncated normal random numbers we will use 

will be from the discrete set of points in (0.00 -.99), so 

/\ 

U 

that really we are only interested in values of o which 

TFT 

are discrete points in (0.00 - O.J4.9). This would imply that 

a maximum of 5 1 graphs would be needed to cover all possible 

cases in which we have an interest. 

Figure III-3 is a set of graphs of the G.P.D. F, for 

R 

truncated normal distributions as described, and for o, the 
mode, equal to 0.00, 0.10, 0.20, 0.30, 0 • i+0, and 0.50. An 
examination of the curves indicates that as the mode approaches 
zero, the successive curves lie closer together. As an ex- 
ample, if we had a uniform random number of .5, and converted 
it using curve #6, we would get a truncated random number 
of .19. If we used graph #5# we would get .22. Now graphs 
#5 and 6 have very near maximum separation of the point used 
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in the example. This means that for an interval that covers 
10 possible positions of the mode, the resolution in random 
number conversion is 3. Apparently, for this interval, we 
do not need graphs for all possible positions of the mode. 

In fact, the graphs for 0.00, 0.05 and 0.10 should be adequate 
for our approximations. 

A similar analysis for successive pairs of curves in- 
dicates that we can get adequate coverage of the cases we 
are interested in by using graphs for 

0.00 

0.05 

0.10 

0.13 

0.16 

0.20 

0.22 

0.24 

0.26 

0.28 and 

0.30 (.01) 0.40, or a total of 

30 graphs. 

If we are only interested in a "fair' 1 approximation, 
then the use of every other curve should suffice, and a total 
of 15 graphs would be adequate. 

In order to construct the necessary graphs, a table of 
the normal integral is required. In the following discussion, 
the term $(a) represents the area under the normal integral 
from - 00 to a J /. e • 
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$(a) ~ J vk- e and 

-CD 

_ t z 

/hf e ^ dt = £(h) - $(<*) 

w 



A property of the normal integral is that 
*(-<») = /- 



The mean ('Y ) of the normal integral must be zero to 
use the tables in customary form. Therefore, we will perform 
our computations in such a manner as to effectively consider 
-*<' = 0. Values less than the mode (which is equal to -y ) 

will therefore be considered as negative. We can accomplish 
this by measuring all our distances in terms of ^ and 
measuring from the point which is the mode. Under our as- 
sumptions, we determine 4 by 

3 . 5 d = /- rr? where AT7 = 

/Z 



or 



a /-/77 

<J. s' 



If our horizontal axis is Y , then all values on 
scaled to some value of <6 , with 0 at m on 

ample, m = .3 on X we would have 



Y can be 
If for ex- 



0 •/ «2 .3 > S & *7 .8 <3 AO vo/ues 

-j 1 1 1 1 1 1 1 1 1 i x 

-/>& -AO - O-S O O S AO AS" 2-0 2 -S So 3 -f I/q/uCS 



84 



















'J 

' 



V 

, • .< 4 * ■« * * » ,m ! •» " — » •**» 

, v :>«* %• *ym m • - * * 

* 

1 

4 

Mi 4 jl* • . * •• * s - 



- 









A 



U- 









S' 






Next, we determine A, the area under the normal density 
from our lowest <4 value (equal to zero on the X scale) to 
our highest <3 value (equal to 1 on the X scale). In the 
example above 

A = J Vzv e = J(3 *r<s ) -$(-/■*£) 

-AS 6 

A - <f(3.s<s\+ £(/.&*)-/ 

Using a table of the normal distribution, 

A = .9330 

For successive points corresponding to 0.00, 

0.05, 0.10, etc., on the X scale, compute using the table, 

I ft 6) - JFC-/-SZ) _ 

A ~ A 

The values obtained are points on the C.P.D. F. curve 
at the corresponding value of X , and can be used to draw 
the desired curve. 
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SECTION $ 



Extensions to the Algorithm 

Throughout the development of the algorithm continued 

efforts have been exerted to develop a tool with sufficient 

flexibility that the user can best fit to the individual 

needs of his organization. Obviously there will be cases 

where it is impractical or impossible to secure professional 

estimates of t^, t and t . 

o' m p 

Ilany war games evaluate weapons systems that are still 
in the process of development. In these cases little more 
than a vendor's estimate of the equipment's performance times 
is available. In other cases new techniques and tactics are 
evaluated in areas where there is no reservoir of experience 
to offer assistance in establishing the appropriate delays 
for performance or reaction times. 

On the surface it xfould seem that in this situation the 
war gammer has little choice but to resort to arbitrary or 
fixed time delays. 

With the advent of ever more sophisticated weapons, the 
ability or inability of an opponent in a war game to use a 
particular piece of equipment in a given situation can con- 
ceivably affect the outcome of the entire game. Under such 
circumstances, the controllers of the game must consider 
every avenue that affords a closer approximation to reality. 

To restate the underlying theme of this thesis; "there 
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exist situations where the use of random time delays is not 
a technical nicety, but a practical necessity to insure va- 
lidity. ” 

How could the algorithm of this thesis be used in situ- 
ations such as those described above. Let us devise an ex- 
ample. 

Suppose that we as war gammers are conducting a game in 
which one of our perimeter defense weapons is an as yet un- 
evaluated mobile rocket launcher. Let us further suppose we 
have only the contractor's estimate of the average time to 
put this weapon into operation; say it is six hours. It 
would Indeed be a rare event if this estimate was accurate. 
But still, some estimate must be accepted as the best 
available at this stage of development. 

One can rather easily decide on some most optimistic 
time short of which, due to physical limitations, the set- 
up of this equipment is virtually impossible. Let us assume 
this time (t 0 ) is four hours. In a similar fashion we can 
establish our most pessimistic time (t p ) at say 12 hours. 
Lastly, we can accept the contractor's estimate of t m as 
six hours or adjust it as our common sense dictates. We now 
have the necessary estimates to produce a random time delay 
of between 4 and 12 hours where, of course, the bulk of our 
generated times will tend to fall near t . 
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APPENDIX 1 



THE AIR BATTLE ANALYZER 

The Air Battle Analyzer was devised by Dr* M. 0. Waddel, 
Applied Physics Laboratory, Silver Spring, Maryland, It has 
recently been used by the United States Marine Corps in a 
study undertaken to evaluate the Air Defense capabilities of 
the Marine Corps Expeditionary Force. 

The Air Battle Analyzer is a scenario type, free play 
war game in which each opponent is given a specific capa- 
bility and plan of action. It attempts to provide, and 
display, a readily accessible means for recording and dis- 
playing chronologically the principle movements and operations 
of the different surface and airborne units involved in a 
battle. This display points up the interactions between 
different units, in particular their ordering in time. 

The physical equipment associated with the Air Battle 
Analyzer consists of a plotting and display chart, several 
tools for measuring and plotting on the displays, and several 
nomographs . 

The chart combines three display plots, a range-azimuth 
plot, a range -altitude plot, and a range-time plot. The 
three plots have a common (horizontal) range scale with 
range origins on a common vertical, to facilitate reference 
from one plot to another. This combination of plots allows 
for the simultaneous visualization of the elements of a 



90 






' 

, 

1 M 














hypothetical air/ground or sea battle in the four dimensions 
of range, bearing, altitude and time. 

Of the various tools and nomographs associated with the 
Air Battle Analyzer only three will be described as being 
germane to this discussion; the Fire Power Analyzer, the 
Range of First Detection nomograph and the Target Speed 
Scaling nomograph. 

FIRS POWER ANALYZER 

A precise estimate of firepower of a surface-to-air 
missile system requires detailed examination of the particu- 
lar system operation. The cycle time between successive 
intercepts, and the intercept envelope depend, often in a 
rather involved manner, upon a large number of system 
parameters. However, a reasonable understanding of missile 
system firepower can be obtained from consideration of a 
hypothetical system. To this end, a simple missile system 
employing one launcher and one director has been defined. 

The Firepower Analyzer is a tool for estimating the 
firepower of a launcher defined as above, employing a 100 n.m. 
missile, and engaging Mach 1.0 aircraft in a wave attack. 

This speed was chosen arbitrarily. When the Analyzer is 
placed on the range-azimuth plot of the display chart with 
center at the missile ship and axis parallel to the attack 
path, the various curves of the Analyzer intersect the 
attack path at successive intercept points. To determine 
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the firepower for any maximum intercept range of 100 n.m. or 
less, first determine the points of first and last intercept 
along the attack path. The first of these will be limited by 
maximum missile range or by range of detection, etc. The 
latter point is determined by t^eapon release point. To esti- 
mate the single launcher system firepower, place the Analyzer 
as indicated above and count the intercepts from 100 mile 
range in to the last intercept point, and then subtract 
those occurring beyond the assumed first intercept point. 

In computing the difference, that is, in counting intercepts, 
it is appropriate to interpolate for fractional intercepts, 
inasmuch as the computation is directed toward estimating 
the expected number of intercepts. To the above difference, 
one more intercept should be added; the one occurring at 
maximum intercept range, to obtain the total firepower. 

FIRST INTERCEPT RANGE 

The range of first intercept by a surface-to-air missile 
system may be limited by target detection range and subsequent 
delay for decision, etc., rather than by the maximum range of 
the missile. In this event the range of first intercept is 
the detection range less the distance traveled by the target 
during the delay and the missile time of flight, and so is 
dependent on target speed. 

This game uses a nomograph, enclosure 1, to facilitate 
the calculation of first intercept range. Enter the right 
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hand abscissa with the detection range, move vertically up- 
ward to the appropriate delay; then move horizontally to the 
left to the proper target and speed; next move vertically 
downward to the same detection range, and read off first 
intercept range from the loiter ordinate. The dashed linos of 
the figure trace the above path for a detection range of 
150 n.m, (A), delay of 3 min. (B), and target speed of Mach 
1,2 (G), giving first intercept range of 7k- n.m. (D and E). 
TARGET SPEED SCALING 

The firepower (expected number of intercepts) of a 
surface-to-air missile battery depends critically upon target 
speed inasmuch as fast targets will pass through the missile 
envelope in shorter time and so be faced with fewer missiles 
than slower targets. To facilitate scaling to various target 
speeds the Air Battle Analyzer uses a nomograph, enclosure 2. 
Enter the right hand abscissa with the number of intercepts 
(A N) obtained from the Firepower Analyzer, move vertically 
upward to the desired target speed; then move horizontally 
to the left, reading off the firepower from the ordinate. 

The dashed lines of the figure trace the path for a differ- 
ence An = 14 (A), and target speed of Mach 0,9 (B), giving 
firepower of 16 (0). The left hand portion of the nomograph 
pemits ready conversion of firepower to an estimate of ex- 
pected number of kills. From the firepower reading on the 
ordinate, continue to the left to the appropriate kill 
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probability, then move vertically downward, reading off the 
expected kills from the left hand abscissa. Again referring 
to the dashed lines of the figure, firepower of 16 (C) and a 
kill probability of 0.6 (D) give expected kills of 9.6 (E). 
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APPENDIX 2 



SAMPLE QUESTIONNAIRE 



You are asked to furnish your best independent, pro- 
fessional estimate of the time it takes to accomplish 
certain assignments. The result as furnished by you 
and other Marine officers will be used in an attempt to 
determine whether such a method can yield more accurate 
estimates than are currently provided by PM's, TM‘s and 
various planning publications. 

Enclosure (1), page 2, lists the events and space to 
fill in your estimates. Should you feel that certain 
events are so specialized in nature that you, in the 
light of personal experience, duty assignments, etc., 
cannot give a reasonably valid estimate, you may omit 
estimates for those events. It is recognized that 
each of these events depends on a wide variety of cir- 
cumstances. I\ 7 o attempt to completely describe the 
situation can be successful. Therefore, it is re- 
quested that you make the best estimate possible in the 
light of your interpretation of each problem posed. The 
time required for the following is offered in explanation 
of the various times you are asked to estimate. 



Most likely time 



The time required to 
accomplish the action 
tinder "normal" cir- 
cumstances. "Normal" 
to be construed as 
those day to day con- 
ditions you would 
expect to encounter 
under circumstances 
such that the described 
action could logically 
be assigned 



Most optimistic time 



The minimum time required 
to accomplish the action 
of all circumstances, 
personnel reaction, weather, 
terrain, etc., combined to 
work in your favor. 



bst pessimistic time 



The time required to ac- 
complish the action under 
extreme adverse conditions. 
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and poor response on the 
part of the participants. 

The conditions are not so 
severe as to preclude ac- 
compli shinent of the mission, 
but weather, fatigue, human 
error, etc., all exercise 
strong degrading effects. 

Note: Unless otherwise indicated, times are to be from 

receipt of assignment until accomplishment, and 
should reflect planning, reconnaissance, and 
coordination requirements. 

Estimates are desired to the nearest 1 5 minutes 
unless otherwise indicated. 

Percentage of times the event will occur within 
h ours/minutes of the most likely time. 

If the event were repeated a large number of times, 
this is your estimate of the per cent that require 
a time within +/- hours/minutes of your estimate 
of the most likely time. 

1st Event 



deserve battalion to organize and occupy combat outpost 
to include preparation of hasty defense and preregistration 
of supporting arms. Assume foot march time from reserve 
position to OP site is 2 hours. 

Most likely time 

Most optimistic time 

Most pessimistic time 

Percentage of times event would occur within 30 minutes of 

the most likely time 

2nd Event 



To mount out "ready" B LT, from Marine Corps Base to afloat 
status . 

Most likely time 

Most optimistic time 
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lost pessimistic ti ie 

Percentage of times event would occur* within 60 minutes of 

the most likely time 

3rd Event 

Infantry company to conduct physical fitness testing using 
current test. Time from first man starts until last man 
finishes. Assume no parallel facilities. 

Most likely time 

Most optimistic time 

Most pessimistic time 

Percentage of times event would occur within 30 minutes of 

the most likely time 

kth 3vont 

To heli-lift rifle company in vicinity of Helipad. Five 
men heliteams. Pad capacity - 3 helicopters. Sufficient 
helicopters airborne in vicinty to move company. Time from 
touch dox*m of first flight until last flight is airborne. 
(Time to nearest 5 minutes) 

Most likely time 

Most optimistic time 

Most pessimistic time 

Percentage of times event would occur within 5 minutes of the 

most likely time 

5th Mvent 



Reserve rifle company to provide squad size daylight 
reconnaissance patrol. Time from receipt of order to departure 
of patrol, but to include briefing, rehearsal, and check. 

Host likely time 

Most optimistic time 
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Most pessimistic time 

Percentage of times event would occur within 1 5 minutes of 
the most likely time 
6th Event 

Time to establish Infantry Battalion OP. Time from arrival 
of first elements at CP site until CP is established and 
functioning adequately. 

Most likely time 

Most optimistic time 

Host pessimistic time 

Percentage of times the event will occur within IS minutes 

of the most likely time 

7 th SVent 

Time to prepare one platoon from an PUP unit to conduct 
formal guard mount, to include as you feel might be required, 
classroom training, rehearsals, and preliminary inspection. 

Most likely time 

Most optimistic time 

Host pessimistic time 

percontage of times the event will occur within 30 minutes 

of the most likely time 

6th Bvent 

Tine to establish one platoon road block in hasty defense, 
800 meters forward of parent company *s MLR position. 

Host likely time 

Most optimistic time 

Host pessimistic time 
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percentage of times the event will occur within 1 $ minutes 

of most likely time 

9th Hvent 

Time to establish field mess. Time from arrival of person 
nel and equipment until mess is set up and ready to feed "A” 
type rations. 



Host likely time 

Most optimistic time 

Host pessimistic time 

Percentage of times the event will occur within If? minutes 

of most likely time 

10th Fvent 



You are a Forward Air Controller directing a live ordnance 
air exercise on Drowns Island. You are also acting as Division 
..'.p . Assuming that you are booked in by radio net to 

a group headquarters at Cherry Point. Hhat is your ostimate 
of the time required from ti^e of radio contact with the a/c 
group headquarters until a flight of four jet a/c are on 
target under your direction as FAC. Assume the group has 
alerted a squadron to have pilots and armed a/c standing by. 
assume a ten minute pilot briefing by the squadron, and conside 
ordnance plug in. (To nearest £ minutes) 

Most likely time _____ 

„;ost optimistic time 

Most pessimistic time 

Percentage of times the event will occur within 15 xuinutes of 
the most likely time 



Your rank 



Specialty 



Years of Service 



